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FORWORD 


This  is  the  final  report  under  Contract  F33615-71-C-1289  for  the 
Aerospace  Research  Laboratories,  Air  Force  System  Command,  United  States 
Air  Force.  It  covers  the  period  from  1  February  1971  to  30  November  1972. 
Two  journal  articles  have  also  been  written  under  the  sponsorship  of  this 
contract  "Electrical  Conductivity  from  Electric  Arc  Measurements",  by 
R.  S.  Devoto  and  D.  Mukherjee,  to  be  published  in  Journal  of  Plasma 
Physics  and  "Measurement  of  Argon  Electrical  Conductivity  at  High  Pressure", 
by  U.  Bauder,  R.  S.  Devoto  and  D.  Mukherjee,  submitted  to  Physics  of 
Fluids.  The  senior  author  of  this  report  is  now  located  at  the  School 
of  Mechanical  Engineering,  Georgia  Institute  of  Technology.  Technical 
Monitors  for  this  report  were  Dr.  U.  Bauder,  Thermodynamics  Laboratory, 
and  Capt.  Allen  M.  Hunter,  II,  Plasma  Physics  Laboratory,  Aerospace 
Research  Laboratory. 
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ABSTRACT 


Equations  are  developed  which  govern  the  behavior  of  a  wall-stabilized 
electric  arc  ii1  which  heat  transfer  takes  place  via  conduction  and 
radiation  emission  and  reabsorption  within  the  heated  gac.  The  differential 
approximation  is  used  to  describe  the  radiative  transport.  Several 
msierlcal  schemes  are  described  which  were  used  in  an  attempt  to  Integrate 
these  equations;  none  were  found  to  work  for  all  values  of  the  absorption 
coefficient.  A  discussion  of  the  development  of  higher-order  differential 
approximations  is  also  Included. 
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INTRODUCTION 


An  important  need  for  application  to  hypersonic  wind  tunnels  is  a 
Bource  of  high-pressure  (to  1000  atm)  high  temperature  (to  15000°K)  air. 
The  high  pressure  electric  arc  would  satisfy  this  need;  its  properties 
are  as  yet  only  imperfectly  known.  For  th^s  reason  a  program  was  under¬ 
taken  to  compute  properties  such  as  temperature  profiles,  power  input, 
etc.  of  a  particular  shape  of  high-press jre  arc,  namely  a  wall-stabilized 
cylindrically  symmetric  arc. 

At  high  pressures,  account  must  be  taken  of  energy  transport  via 
radiation  reabsorption  within  the  arc  as  well  as  by  the  traditional 
thermal  conduction.  However,  the  exact  theory  of  radiative  transport  is 
complicated,  and  therefore  difficult  to  apply  in  real  problems.  For 
this  reason  many  people  have  employed  the  so-called  differential 
approximation.  The  computation  of  arc  properties  with  this  approximation 
is  described  in  Section  II.  Because  of  doubts  about  the  accuracy  of 
this  approximation  for  certain  values  of  the  optical  depth,  we  have  also 
Investigated  higher-order  methods.  This  work  is  described  in  Section  III. 
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SECTION  II 

LOtfEST-ORDER  DIFFERENTIAL  APPROXIMATION 


For  a  cylindrical ly  symmetric  arc,  the  gas  energy  equation  and  the  D1 
equations  for  radiative  transfer  can  be  written  In  the  form 


l 

r 


(1) 


dT 

dr 


f  (qT-0R) 


(2) 


J. 
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,  * 
cx  (u 

V  V 


U  ) 
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(3) 


(4) 


T  R 

where  q  and  q  are  the  total  and  radiative  contributions  to  the  heat  flux, 

x  is  the  effective  absorption  coefficient,  u  is  the  radiation  energy 
v  *  v 

density,  and  u^  is  its  equilibrium  value,  given  by 


u*  .  2ri>sil _ i _ 

V  c3  .WkT.j 


(5) 


Other  symbols  have  their  usual  meaning. 

Boundary  conditions  for  these  equations  can  be  written  in  several 
equivalent  ways.  In  the  present  case  we  take 


T(0) 

-ta 

(6) 

qT(0) 

-  0 

(7) 

qR(0) 

-  0 

(8) 

T(R) 

(9) 

CU  (R) 

-  2qR  (R) 

(10) 

v  v 
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Equations  (6)  -  (9)  are  obvious  on  physical  grounds;  Eq.  (10)  arises  from 
the  necessity  for  continuity  of  heat  flux  at  the  arc  outer  boundary. 

He  see  that  these  are  lour  coupled  first  order  differential  equations, 

but  with  five  boundary  conditions,  seemingly  one  too  many.  One  of  the  boundary 
conditions  is  used  In  determining  the  value  of  E  ,  so  this  Is  an 
eigenvalue  problem.  In  an  alternative  formulation,  E  Is  given  and  T(0) 

Is  computed  as  part  of  the  solution  ,so  Eq.  (6)  Is  no  longer  needed. 

The  above  equations  can  be  transformed  Into  a  form  more  convenient 

for  solution  if  we  define  new  dimensionless  variables 
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Equations  (1)  -  (4),  (6)  -  (10)  then  become 
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u(l)  -  T^/Ta 

(25) 

4z  (D 

w  (1)  -  -r^ - 

v  3t.  * 

(26) 

at  x  ■  0  ,  the  limiting  forms  of  Eqs.  (18)  and 

(20) 

are 

du  j  \ 

/* 00 

dv  c  1 

UuA 

dz 

dx  "  ’  \ 

r '  l 

\A 

dx  d 
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(18a) 


(20a) 


which  must  be  used  near  x  =  0  to  avoid  the  singularities  in  Eqs.  (18) 
and  (20). 

Preparatory  to  solving  the  full  set  of  equations  above  we  can  look  at 
several  limiting  cases. 

1.  Negligible  Radiation  -  In  this  case  the  equations  become  the  well-known 
Elenbaas-Heller  equation.  In  the  form  of  the  equations  used  here,  the 
mathematical  problem  is  a  two-point  boundary  value  probloa  l.e.,  the  two 
first-order  (or  one  second-order)  equations  must  be  solved  subject  to  condi¬ 
tions  at  two  boundaries.  In  the  present  case  we  can  numerically  Integrate  the 
differential  equeMons  (18) -(19)  as  an  initial-value  problem  using  the 
boundary  conditions  (22) -(23)  as  initial  values.  A  value  is  assumed  for 
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o  (s^y  crj )  and  the  equations  Integrated  to  x  •  1  .  If  £q.  (25)  is  satis¬ 
fied  to  within  a  small  constant  t  ,  then  or^  is  the  eigenvalue  and  the  com¬ 
puted  solution  is  the  desired  one.  In  general,  Eq.  (25)  will  not  be  satisfied 
and  wf  must  guess  a  new  value  for  a  >  say  05  •  The  procedure  for  computing 
the  new  value  of  a  must  be  chosen  with  care  since  a  great  deal  of  computation 
time  can  be  used  while  integrating  the  differential  equation  from  x  *  0  to 
x  «■  1  with  new  values  of  Q  .  In  the  present  case,  after  the  initial 
guess  ,  e»2  was  computed  from 

!0.01  #1  ,  u ( 1 ; Qfl ) c  T^/Ta 

(27) 

-0.01  ,  u ( 1 ; >  1^/T 


and  further  cr  from 


Vi 


or  i  .  a  . 

1  ~  V 


u(l ; or  ,)  -  «(l;o  ,) 

V-l  \Tl 


u(l,a  )-  VTa 


(28) 


The  dependence  of  u  on  er  1»  explicitly  indicated  above.  An  example  of 

v  4o 

a  temperature  profile  confuted  for  argon  at  1  atm  pressure  with  T  -  10  K 

3 

and  ■  10  °K  is  given  in  Fig.  1.  Shown  is 


T  -  \ 


vs.  x  rather  than  u(x)  . 

2.  Cons tant  Effective  Absorption  Coefficient  -  In  this  case  Eqs.  (18)  -  (22) 
reduce  to  the  following 


du 

dx 


*) 


(29) 


P-  -  0  X-  (30) 

dx  afl 
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dw 

dx 


z 

X 


(31) 


dz 

dx 


(w* -w) 


(32) 


with  correspondingly  simplified  boundary  conditions.  N  Is  the  so-called 
conduction -radiation  parameter  and  is  given  by 


N  « 


(33) 


where  d  Is  the  Stefan-Boltzmann  constant,  and  Is,  crudely  speaking,  a 
measure  of  the  relative  importance  of  conduction  and  radiation  heat 
transfer  in  the  arc.  Other  quantities  are  r  =  Rh  (a  constant)  and 
w*  ■  u\ 

As  a  preliminary  to  Integrating  the  arc  equation  with  varying  k  . 
the  solution  of  these  equations  was  programmed  for  the  IBM  360.  Important 
to  the  program  are  two  major  parts:  (1)  the  method  for  numerically  in¬ 
tegrating  the  differential  equations,  (2)  the  method  for  satisfying  simul¬ 
taneously  the  boundary  conditions  at  x*0  and  x  -  1. 

In  view  of  the  fact  that  the  differential  equations  probably  become 
stiff  for  T  ;>  1  ,  it  is  necessary  to  use  a  technique  applicable  to  such 
equations.  A  program  for  integrating  systems  of  stiff  or  non-stiff  equations 
has  been  recently  published  by  Gear1  and  was  used  for  the  present  case.  It 
was  found,  however,  that  under  certain  conditions  the  Gear  program  would  not 
integrate  the  equations  and  exited  with  an  error  message.  Several  months  of 
effort  failed  to  correct  this  problem  and  it  was  only  recently  learned  that 
this  program  has  been  rewritten  by  Gear  to  remove  this  error.  The  new  Gear 
program  has  since  been  incorporated  into  the  program  and  has  given  no  more 
difficulties. 

The  second  problem  listed  above  has  proved  much  more  difficult  and  was 
not  resolved  in  this  work.  For  Eqs.  (29)-(32)  there  are  boundary  conditions 
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to  be  satisfied  by  u(as  in  the  Elenbaas -Heller  equation)  and  by  z  and  v  . 

To  integrate  the  equations  as  an  initial-value  problem  we  must  guess  a  *nd 

w(0),  integrate  to  x  -  1  and  see  if  Eqs.  (25)  and  (26)  are  satisfied.  If 

they  are  not  satisfied,  then  a  and  w(0)  must  somehow  be  adjusted  in  an 

iterative  process  to  do  so.  As  a  first  attempt  to  iterate,  after  initial 

guesses,  a  and  v(0)  were  adjusted  as  in  Eqs.  (27) -(28).  Tins  procedure 

was  successful  for  small  t  ,  but  did  not  converge  for  t>0.01.  Profiles 

of  9  vs  x  for  t  -  0.0025,  0.01  are  shown  in  Fig.  1  for  an  argon  arc  at 

IA  ■  10  °K  ,  p  ■  1  atm.  We  note  that  the  temperature  profile  is  depressed 

near  the  axis  due  to  radiative  losses,  as  we  would  expect. 

Since  the  above  procedure  was  not  successful  for  t>  0.01  ,  a  new 

method  was  sought.  It  turns  out  that  the  radiation  density  at  the  arc 

center  can  be  evaluated  from  an  exact,  fairly  simple  two-fold  quadrature 

once  the  temperature  profile  (or  an  approximation  thereto)  is  known.  In  the 

2  3 

notation  of  this  report  the  general  expression  reads  ' 


«  (0) 


x  (x)  w*(x)  D.  T  (0,x) 

V  V  1  V 


dx 


(34) 


where 

t  <0,x)  - 
v 

is  the  optical  depth  at  frequency  v  and  position  x  and  D^(x)  is  an 
Integrated  Bessel  function, 


/ 


X  (x')dx' 

V 


(35) 


Kc(y)dy 


(36) 


For  constant  x  Eq.  (34)  reduces  to 


v(0) 

with  t  *RX 


u*(x)  Dj(Tx)dx 


(37) 
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The  new  interation  procedure  is  as  follows: 

(1)  Compute  u(x)  without  radiation  or_  feed  in  assumed  u(x); 

(2)  Compute  w(0)  from  Eq.  (37); 

(3)  With  '’(0)  from  s*ep  (2)  and  guessed  or^  integrate  arc  equations 
until  convergence  is  obtained  on  a  ,  i.e.,  the  eigenvalue  a  for  the 
w(0)  computed  in  (2)  is  obtained; 

(4)  Recompute  w(0)  with  Eq.  (37)  and  repeat  (3); 

(5)  After  w(0)  no  longer  changes,  check  to  see  if  Eq.  (26)  Is  satisfied. 
If  so,  this  is  the  desired  solution;  if  not,  then  apply  a  scheme  similar  to 
Eqs.  (27) -(28)  to  both  or  and  w(0)  until  convergence  is  obtained.  Steps 
(l)-(4)  should  ensure  that  the  initial  or  and  w(0)  are  close  enough  to 
the  true  solution  to  obtain  convergence. 

This  procedure  was  tried  for  t  >  0.01  but  it  also  did  not  converge. 
The  reason  for  the  failure  probably  lies  in  the  value  of  w(0)  computed 
from  Eq.  (37).  This  is  tin  exact  value  if  u(x)  is  exact.  However,  the 
solution  is  an  approximate  solution  which  will  presumably  deviate  from 
the  exact  solution.  We  thus  expect  that  w(0)  computed  from  Eq.  (37), 
even  with  the  exact  u(x),  will  not  be  the  correct  w(0)  for  this 
problem.  Note  furthermore  that  a  small  error  in  u(x)  is  multiplied 
fourfold  in  the  computation  of  Eq.  (37). 


SECTION  III 


HIGHER-ORDER  DIFFERENTIAL  APPROXIMATIONS 


Considerable  effort  has  been  directed  towards  obtaining  solutions  of 
the  radiative  transfer  equation  and  its  close  counterpart  in  neutron 
transport  theory.  A  thorough  literature  search  carried  out  indicates 
that  this  effort  can  be  subdivided  into  the  following  major  categories: 

a.  Exact  solutions 

b.  Purely  numerical  solutions,  including  the  Monte  Carlo 
t  echr ique 

c.  Approximation  schemes 

Because  of  its  complicated  integro-dif ferential  nature,  exact  solutions 
have  been  obtained  only  in  very  simple  cases  with  restrictions  on  the  absorption 
coefficient  and  geometry  of  the  system.  Together  with  the  second  category 
(numerical  solution)  the  exact  solutions  are  of  interest  in  so  far  as  they 
provide  a  check  on  the  results  from  approximate  techniques  of  solution. 

Most  of  the  approximation  schemes  are  related  to  one  or  the  other  of 

two  general  methods  for  approximating  the  angular  dependence  of  the  radiation 

intensity  —  the  spherical  harmonlcb  method  (or  the  method)  and  the  discrete 

ordinates  method  (or  the  S  method).  Interest  in  the  latter  method  has  been 

n 

very  limited  and  ve  have  primarily  concerned  ourselves  with  the  and 
related  methods  in  our  effort. 

The  straightforward  spherical  harmonics  approximation  consists  of  ex¬ 
panding  the  radiation  intensity  I  in  terms  of  spherical  harmonics  Y™ 

(which  reduce  to  Legendre  polynomials  in  planar  and  spherically  symmetric 
geometries  and  to  associated  Legendre  polynomials  in  a  cylindrically  symmetric 
case),  substituting  into  the  transfer  equation  for  I  ,  multiplying  the 
resulting  equation  by  Y®  ,  the  conjugate  spherical  harmonics,  and  integrating 
over  the  conq>lete  solid  angle  range.  Use  of  the  orthogonality  and  recurrence 
relations  for  the  spherical  harmonics  leads  to  an  infinite  series  of 
ordinary  differential  equations  in  terms  of  the  coefficients  of  the  spherical 
harmonics  expansion,  ,  which  are  functions  of  position  only  and  are  thereby 

related  to  the  quantities  of  interest,  viz.,  the  radiation  intensity  and  heat 
flux  vector. 
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Hie  P  approximation  consists  of  setting  the  coefficients  A  ” 
n  1 

for  ni  N+l,  thus  leading  to  a  finite,  determinate  set  of  ordinary  differ¬ 
ential  equations.  Hie  first  approximation  or  Pj  approximation  is  often 
referred  to  as  the  differential  approximation. 

Although,  as  formulated,  this  theory  has  no  explicit  restrictions  on 
either  absorption  coefficient  or  geometry,  it  turns  out  that,  as  a  practical 
matter,  the  P^  or  P^  approximations  fail  whenever  the  directional 
distribution  of  the  radiation  intensity  becomes  singular  or  when  any  volume 

element  in  the  ga3  is  not  traversed  by  photons  coming  from  all  directions. 

Such  is  the  case  for  the  radiation  field  in  curvilinear  geometries  such 

as  concentric  spheres  or  cylinders. 

Several  techniques  for  overcoming  this  problem  have  been  published  in 

U 

recent  years.  One  of  the  more  promising  of  these  is  that  by  Traugott  in 
which  a  modified  spherical  harmonic  technique  is  applied  to  the  radiation  field 
between  concentric  spheres.  Instead  of  developing  the  straight-forward 
spherical  harmonics  approximation,  Traugott  develops  a  finite  set  of  moment 
equations  where  the  moments  are  those  of  I  with  respect  to  direction  cosines 
of  0  .  As  a  system  of  moment  equations  is  indeterminate,  it  is  necessary 

to  obtain  closure  relations  between  the  moments  with  the  aid  of  a  spherical 
harmonics  expansion  of  the  intensity.  It  is  in  these  closure  relations  that 
the  modification  is  made  by  requiring  that  these  relations  hold  not  only  for 
the  case  of  isotropy  but  also  for  a  unidirectional  (radial)  intensity  field. 

The  attractiveness  of  this  approach  lies  in  the  fact  that  it  is  possible 
to  incorporate  into  the  approximation  scheme  physical  arguments  relevant  to 
the  particular  geometry  under  study.  With  this  in  view,  our  main  effort  this 
quarter  has  been  towards  developing  general  moment  equations  and  closure 
relations  between  moments  in  cylindrical  geometry.  These  results  are 
summarized  in  the  following  sections. 

1.  Transfer  Equation 

The  radiative  transfer  equation  for  a  non-scattering,  steady-state  gas 
in  local  thermodynamic  equilibrium  can  be  written  as  ^ 

•  $  Iv(?.n)  -  •  Xv(T)  (Iv(r,fi)  -  Bv(T)  )  (38) 

where  (t  is  a  unit  vector  in  the  direction  of  photon  propagation,  1^ 
is  the  spectral  radiation  intensity,  is  the  P  lanck  function  for 
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equilibrium  radiation,  is  the  effective  absorption  coefficient  and  T 

is  temperature.  The  Planck  function  B  is  defined  bv 

v 


B  -  2hv3  1 

'  c2  (exp(hv/kT)-l) 


(39) 


Furthermore,  the  radiative  heat  flux  vector  -  the  quantity  of  main 
interest  to  us  —  is  defined  by 


<1 

v 


i  n  d  n 

V 


(40) 


For  the  sake  of  convenience,  we  shall  suppress  the  functional  dependence  of 
these  quantities  as  well  as  the  subscript  "v"  . 

It  is  convenient  to  rewrite  the  transfer  equation  in  terms  cf 
cylindrical  polar  coordinates  at  the  outset.  Using  the  coordinate  system 
shown  in  Fig.  2  and  assuming  cylindrical  symmetry,  we  find  for  Etjn.(3fi): 


sin©  cos  y  £  -  &  -  -  X(I-B)  (41a) 

Alternatively,  this  can  be  written  in  terms  of  direction  cosines  as 

i  hL-  al  «  .  x(I-B)  (41b) 

r  *r  r  ay 

where  the  direction  cosines  are  given  by 

2  1/2  //o  ^ 

l  «*  sin©  cosy  ■  (1-^  )  cosy  (42a) 

l  -  sin©  sinX  «  (l-^i2)1^2  ain  y  (42b) 

© 

-  COS  0  -  |A  (42c) 

2,  Moment  Equations 

We  define  a  sequence  of  moments  of  the  intensity  as  follows: 

I<)00 '  'C  “n  (43‘) 

■too  -  ^Itrdn  <43b> 

■oio  -  ■yn  <43c> 
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m  nf'iHH!'  r 


i 


Figure  2.  Coordinate  System  In  Cylindrical  Geometry  (c  corresponds  to 
axis  of  cylinder). 
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001 

=  f  n  dn 

</  z 

4rr 

(43d) 

110 

*=  J  \l  l  dn 

4tt  r  * 

(43e) 

101 

=  /  in  dn 

4tt  r  Z 

(43f) 

Oil 

»  f  ll  l  dQ 

4rr  ®  z 

(43g) 

Here,  it  is  understood  that  the  first,  second  and  third  subscripts  refer  to 
the  order  of  l  ,  i  and  i,  »  respectively,  in  the  integrand  on  the  R.H.S. 

r  cp  2 

It  is  easily  shown  that  the  zeroth  order  moment  IQ00  is  proportional 
to  the  radiation  energy  density  u(r)  ,  the  three  first-order  moments 
(I^qq,  IQ10,  Iqqj)  comprise  the  three  components  of  the  radiation  heat  flux 

vector  q  and  the  nine  second-order  moments  are  proportional  to  the  radiation 

pressure  tensor.  Moments  higher  than  second-order  have  no  physical  significance. 

Differential  equations  for  the  moments  are  derived  as  follows:  To  obtain 

the  nth  order  moment  equation  we  multiply  the  transfer  equation  (41b)  by 

f1?  (where  l.  can  be  l  ,  l  or  t)  and  integrate  the  resulting  equation 
i  i  r  cp  z 

with  respect  to  direction  over  the  entire  solid  angle  range.  Thus  the  zeroth 
order  moment  equation  is 


dl 


100 


dr 


—  x 
r  100 


+  x  I 


000 


4  irKB 


(44) 


The  first-order  moment  equations  consist  of  a  set  of  three  equations  as  follows: 


dl 


~  <*200  "  ^20)  +  xl100 


-  0 


(45a) 


dI110  2I110 
- iiii  .  - LiH  +  V  T  .  0 

dr  r  R1010 


dI101  .  *101  .  „ 

~~dr~  +  —+  Kl001  0 


(45b) 

(45c) 


By  continuing  this  process  we  can  generate  an  infinite  number  of  moment 
equations  which,  taken  together,  are  con^letely  equivalent  to  th»  original 
transfer  equation.  An  approximation  consists  of  truncating  this  series  and 
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considering  only  a  finite  number  of  moments  equations.  In  this  case,  however, 
the  n  •mber  of  unknowns  usually  exceeds  the  number  of  equations(and  it  is 
necessary  to  obtain  additional  relations  between  the  moments  so  that  the 
system  becomes  determinate. 

Ue  can  obtain  these  "closure"  relations  as  follows.  Expand  the  inten¬ 
sity  in  a  series  of  associated  Legendre  polynomials: 


I(r,|i,x) 

X-0 


2At-l 

4tt 


P°<u)A°<r) 


J l 

hfI 


( l- m)  1 
( ft-  m)  • 


P1?  (^)cos  V  A“  (r) 


m 


tp 


(46) 


For  the  Nth  order  moment  approximation  scheme  (which  is  obtained  by 
integration  w.r.t.  Q  of  successive  products  of  (41b)  and  j”  ,  m  ■  0,1,2...N) 
we  truncate  the  expansion  (46)  at  l  *  N  ,  by  setting  A™  *  0  for 
/  >  N  +  1  .  Substitution  of  the  truncated  spherical  harmonics  series  for  I 
Into  the  defining  equations  for  the  moments  (43)  and  elimination  of  the 
coefficients  A™  from  the  resulting  equations  give  us  the  necessary  num- 
ber  of  closure  relations  between  the  moments  to  render  the  problem 
determinate.  This  resulting  problem  —  the  moment  equations  plus  the 
closure  relations  —  is  identical  to  the  determinate  system  of  differential 
equations  in  terms  of  the  coefficients  AX  one  obtains  by  carrying  out  the 
straightforward  spherical  harmonics  expansion  (described  in  the  introduc¬ 
tion  to  the  section)  to  order  Pu  . 

For  the  first-order  approximation  scheme,  in  which  equations  (44)  and 
(45)  are  the  moment  differential  equations,  the  closure  relations  are  found 
to  be 


X010  "  X110  *  X101  "  ° 


X200  “  X020  ”  X002  "  1/3  X000 


(47a) 

(47b) 


Substitution  of  (47)  into  (44),  (45)  leads  to 


dl 


4.  i  j  4.  Hi  »  4tt  y  8 
dr  r  l100  +  Hl000  *B 


dl 


000  .  ,  T 
dr" +  3kXioo  0 


(48a) 

(48b) 
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Replacing  1^^  by  ,  the  radial  component  of  the  radiation  heat  flux 
vector  ,  ai  d  Iqqq  by  cu  (where  u  is  the  radiation  energy  density),  we 
reproduce  equations  (3) -(A)  of  the  earlier  section.  In  other  words,  the 
first-order  moment  approximation  scheme  is  identical  to  the  differential 
or  approximation. 

We  have  evaluated  the  differential  equations  and  the  c’osure  relations 
for  the  next  two  higher  moment  approximation  schemes.  These  correspond  to 
the  and  spherical  harmonics  approximations.  As  noted  in  the 

Introduction,  we  are  attempting  to  modify  the  closure  relations  In  these 
higher  approximations  by  incorporating  physical  arguments  relevant  to  a 
curvilinear  geometry.  For  the  P^  approximation,  which  is  being  developed 
at  present,  we  are  modifying  the  closure  relations  so  that  they  satisfy  the 
relationship  between  moments  that  result  when  the  radiation  field  consists 
of  an  undirectional  radial  beam.  In  this  way  the  resultant  differential  equations 
should  be  able  to  more  accurately  predict  the  radiative  heat  flux. 

In  these  approximate  schemes  for  solution  of  the  transfer  equation 
it  is  necessary  to  replace  the  exact  boundary  conditions  by  those  that 
are  consistent  with  the  approximation  being  made.  Two  general  techniques 

have  been  developed  for  this  purpose  and  have  been  applied  extensively  in 

4 

neutron  transport  problems.  One  of  our  primary  goals  was  to  develop  a 
consistent  set  of  boundary  conditions  for  the  higher -order  approximation 
schemes  discussed,  and  then  to  apply  the  differential  equations  plus 
boundary  conditions  to  the  arc  problem. 
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